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The parking-lot model provides a qualitative description of the main features of the 
phenomenology of granular compaction. We derive here approximate kinetic equations 
for this model, equations that are based on a 2— parameter generalization of the statistical- 
mechanical formalism first proposed by Edwards and coworkers. We show that history- 
dependent effects, such as memory and Kovacs effects, are captured by this approach. 

1. Introduction 

The term of "glassy-dynamics" is now commonly used to describe out-of-equilibrium 
systems that display such generic features as very slow kinetics that prevent the system 
from reaching equilibrium in any reasonable experimental timescale, aging phenomena, 
history-dependent processes like hysteresis and memory effects. Among such systems are 
the "not-too-strongly" vibrated granular materials El S I EJ il E| . In the recent years, there 
has been a surge of research activity in this field jy, |7|, |8|, partly driven by the goal of 
providing a statistical-mechanical description of these out-of-equilibrium situations]^, 11(1 

3 E3 B Q S H E3, H Q • 

In this note, we consider an approximate statistical-mechanical description of the parking- 
lot-model (PLM) for vibrated granular materials l2l. la. l20l l2ll l22l l23l . |24| that is based on 
the formalism proposed by Edwards and coworkers [li. lid. Till ll 2t| . Despite its simplicity, 
the one-dimensional model of random adsorption-desorption of hard particles (PLM) has 
the merit of being a microscopic, off-lattice model that mimics many features of the com- 
paction of a vibrated column of grains. It also has, we hope, a didactic value as to the 
nature of several canonical characteristics of "glassy dynamics" . Many of the properties 
of the model, that can be obtained either a naly tically or by computer simulation, have 
been already described in the literature 0, H HQ EBB SB. Our main focus 



here is on memory effects, inclu ding the so-called Kovacs effect first observed in glassy 
polymers j^ - " 



see also[29,, |30j, |31|, |32|), and on the ingredients that are needed in a 



statistical-mechanical description to account for such effects. 



2. The model and its properties 



The parking-lot model is one-dimensional process in which hard rods of length a are 
deposited at random positions on a line at rate k + and are inserted successfully only 
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if they do not overlap with previously adsorbed particles; otherwise they are rejected. 
Moreover, all deposited particles can desorb, i.e., be removed from the line at random 
with a rate k_. For convenience, the unit time is set to l/k + , and the unit length to 
a. With this choice of units, the only control parameter in the model is K = k + /k_. 
When desorption is forbidden (k_ = 0), the model corresponds to the purely irreversible 



one-dimensional random sequential adsorption (RSA) process [33, 34], also known as the 
car parking problem, and all the properties of the system as a function of time can be 
obtained exactly. Connection to the compaction of a vibrated column of grains is made 
by considering the line as average layer (a 2-dimensional model would be more realistic, 
but the qualitative behavior would not be altered), the time as the number of taps, and 
1/K as the tapping strength that controls the fraction of particles ejected from the layer 
at each tap. 

When 1/K is not strictly equal to zero, adsorption and desorption are competing mech- 
anisms that drive the system to a steady state corresponding to an equilibrium fluid of 
hard rods at a constant activity 1/K. All the properties of the steady-state can also be 
obtained exactly. 

The compaction kinetics of the parking-lot model at constant K is described by 



dp 
dt 



= *w-f (i) 

K K 



where p(t) the density of hard rods on the line at time t and $(£) is the fraction of the 
line that is available at time t for inserting a new particle, i.e., the probability associated 
with finding an interval free of particles (a "gap") of length at least 1. The quantities p 
and $ can be obtained from the 1— gap distribution function G(h,t) which is the density 
of gaps of length h at time t via a number of "sum rules" : 



Pit) 
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/ dhG(h,t) = l- dhhG(h,t), (2) 
Jo Jo 



oo 



$(t) = / dh(h-l)G(h,t). (3) 



The 1— gap distribution function G(h, t) obeys a kinetic equation that also involves 2— gap 
distribution function. Similarly, the kinetic equation of the 2-gap distribution involves the 
3— gap distribution function, and so on[25|. The resulting infinite hierarchy of coupled 
equations cannot in general be solved analytically (exceptions are the RSA, when /c_ = 
and equilibrium when t — > +oo). 

The description of the kinetics of the parking lot model at large, constant K (i.e. 
small, constant tapping strength) can be summarized as follows: the density increases 
monotonically during the process, and the kinetics can be considered as a succession of 
four different regimes (see Fig. |TJ): during a first stage, the density increases rapidly until 
a value of around 0.65, and it is followed by an algebraic approach toward the saturation 
density, pjl = 0.747... of the model without desorption (RSA); around this density 
there is a crossover to a still slower l/ln(t) compaction regime that is reminiscent of 
what is exactly observed when 1/K — > 0, and finally there is an exponential approach 
to the steady-state (equilibrium) density with a rate T ~ $ ■ For the two first 
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Figure 1. Logarithmic-linear plot of the density versus time for the parking-lot model at 
a constant tapping strength 1/K = 0.0002 



regimes, desorption has a negligible effect and only the two last regimes are relevant 
for the compaction kinetics of vibrated granular materials (although the approach to 
equilibrium may be prohibited when K is very large). 

Experimentally, it has been observed that granular compaction exhibits history-dependent 
phenomena. In particular, Josserand et. al[A\ have shown that memory effects are seen 
when changing the tapping strength during the compaction kinetics. Similar effects have 
been found in the parking-lot model, when the tapping strength 1/K is switched at a 
given time to a larger (or smaller) value: see Fig. |2] 

The obvious lesson that one can draw such memory effects is that whereas equilibrium 
is fully described by one thermodynamic parameter, the density p, out-of-equilibrium 
situations require at least one additional "thermodynamic" parameter. This parameter 
cannot simply be the tapping strength since the system can be found in states character- 
ized by the same tapping strength 1/K and the same density p, that nonetheless evolve 
in different ways under further tapping with the same strength 1/K : this is illustrated in 
Fig. |2b- A natural candidate for an additional thermodynamic parameter is the available 
line fraction $. One can indeed check that the two states described above do correspond 
to distinct values of $. 



3. Statistical-mechanical formalism with two thermodynamic parameters 

Following the ideas put forward by Edwards and coworkers 0, 11 1, we consider a 
statistical-mechanical description of the system in which all possible microstates character- 
ized by a (small) number of fixed macroscopic quantities are assumed to be equiprobable 
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Figure 2. Memory effect: (a) density versus time when the tapping stength is switched 
from 5.e~ 4 to 2.e~ 3 (lower full curve) at t = 1000 and from 2.e~ 3 to 5.e -4 (upper full 
curve). The dashed curves correspond to a constant strength, 5.e~ 4 (lower curve) and 
2.e -3 (upper curve); (b) Zoom up on the region around t = 1000 when one switches 
from 1/K = 5.e~ 4 to 2.e~ 3 . The arrows denote two points at the same density and same 
tapping strength 1/K. 



("flat" or microcanonical distribution). In addition to fixing the density p, the parameter 
originally selected in the compactivity-based description of granular media by Edwards 
and coworkers, we also constrain the microcanonical ensemble by fixing the available line 
fraction $ to account for the above discussion. 

Denoting by A the total length available for insertion of a particle center (A = $L), 
the configurational integral with the constraints of fixed A, fixed system size L, and fixed 
number of particles N is obtained asfiif 



Z(L, N,A)= £ ... £j[ dh t 6 (l - N - (a - ^ d{K - 1) % - 1) j , (4) 



which can be rewritten as 



Z(L,N,A) 



dz / dyexp ( L 
c J a 



z(l - p) +y® + pin 



z + y{l — exp(— z)) 



(5) 



z(z + y) 

where C and C denote two closed contours. In the macroscopic limit, iV — > oo, L —>■ oo 
A — > oo with p and $ fixed, one can use a saddle-point method to evaluate the integrals 
which leads to 



Z(L,N,A) ~exp(Ls(p, $)) 



(6) 
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p J z z + y z + y(l — e 
$ 1 l-e 



where the entropy density s(p, <3>) is given by 

s(p , «) = (i - P ), + y + P to ( » + v{1 z -^ ( - z)) ) , (7) 

with 2 = z(p, $) and y = y(p, $) solutions of the two coupled equations 

i-* \JL_^i±^, (8) 

(9) 

p z + y z + y(l — e z ) 

The gap distribution functions can also derived along the same lines j^, which leads to 

z(z + y) _ zh , , 

p _ e ztl for h < 1, 

G Ed {h-p) = { z + y, [l z > (10) 
p Z{Z / V) e -^ + y(h-i)) for/l>L 
z + y (1 — e z ) 

It can also be shown that the multi-gap distribution functions satisfy a factorization prop- 
erty, e.g., GEd(h,h f ; p, $) = GEd(h; p,§)GEd(h'; p,Q). Note that the 1-gap distribution 
function is a piecewise continuous function that obeys the exact sum rules, Eqs. (J2J) and 
©. 

A detailed comparison between the predictions of this statistical-mechanical treatment 
and simulation data at the same p and $ can be found in Ref . jlH . The conclusion is that 
although not exact, and even missing some qualitative features in the limiting case of a 
purely irreversible RSA process, the approach provides an overall good description of the 
data. (It becomes of course exact in the steady state since this latter corresponds to an 
equilibrium situation.) 

4. Approximate description of the compaction kinetics 

A quasi-thermodynamic approach is useful because it allows to predict the structure of 
the system (correlation functions) as well as the fluctuations at any given state point (here 
assumed to characterized by p and $). In the cases where phase transitions occur, it can 
also provide interesting constraints and relations between the parameters characterizing 
the phases at the transition or help to determine the limit of stability of the phases in 
mean-field like treatments. However, one still faces the problem of predicting the state of 
the system for a given protocol and a given history: in the simplest case, for a given time 
t and a given tapping strength 1/K. 

We thus want to push the Edwards formalism one step further. Predicting the trajectory 
made by the system in the (p, $) diagram for a given protocol amounts to determine a 
system of equations relating p and $ with t and K. Eq. (0) is one such equation. To 
obtain another equation for the evolution of $, we start for the exact kinetic equation for 
the 1— gap distribution function for h > 1: 



dG(h,t) 



dt 



(h-l)G(h,t)+2 I " dh'G(h',t)- 2G<(h,t K — [ 1 dh'G(h',h-l-h',t), 
Jh+i K K J 



K 
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where G(h, h', t) is the 2-gap distribution function that satisfies the "sum rule" 

dh'G{h,h',t) = p(t)G(h,t). (12) 



f 

Jo 



Note that the approximate 2— gap distribution function GEd{h, h!) satisfies the above sum 
rule. By multiplying Eq. (fTT|) by (h — 1), integrating over h from 1 to oo, and using the 
sum rules, Eqs. (j2J), © and (|T2"|). one obtains an exact equation for the evolution of $: 



dt 



K 



-(1 -p(t) -*(*)) - 



dh(h-l) 2 G(h,t)+ / dh(h-2) 2 G{h,t). (13) 



If one now inserts the approximate expression of the 1-gap distribution function, Eq. ()10|). 
in the above equation, one gets 



<9$ 
~dt 



2(1 -p(t) -$(*)) 
K 



2$(t)- 



_ e -(j/(*)+2(t)) 

y(*) + *(*) 



(14) 



Eqs. (JTJ), (jHJ), and form a closed set of equations whose solution for given initial 
conditions completely characterizes the system and its evolution. 

We first test the accuracy of the above approximate kinetic description in two limiting 
cases for which the exact solution is known: the purely irreversible RSA case {1/K = 0) 
and the approach to the steady state (t — > +oo) for a given rate K. The steady state 
corresponds to the solution = . p °° , i/oo = 0, where is the equilibrium density 

1 poo 

of hard rods at constant activity 1/K: p^ = with $oo = (1 — p^) exp ( ~ p °° ); it 

y 1 Poo J 

is thus the exact result. The approach to the steady state is obtained by linearizing the 
kinetic equations around the equilibrium solution. One then obtains the following coupled 
linear differential equations : 



dy(t) 



dt 

dz{t) 



-2 



(zoo + l)(e 



-2z~ 



-i) + 2(i + 4,K 



dt 



K 



K 



6e z 



Zoo (2c 
2zl 



6 - 8Zr 



2z x 
K2 



-y(t) 



(15) 



Zoo (2<3c 



2z n 



-y(t) - 



Zoo + 1 



(z(t) - Z^). 



(16) 



The associated eigenvalues are all negative and that corresponding to the inverse of 
relaxation time goes for large K as 



3K 



(17) 



The relaxation time is of the same order of magnitude as that obtained in the simple 
adiabatic approximation^ and is much smaller than the exact result (see above). The 
error comes from the inability of the Edwards approximation to account for the non- 
exponential behavior of the 1— gap distribution function for < h < 1. 
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Figure 3. Density versus time for the model without desorption (RSA): the full curve 
corresponds to the approximation, the dashed curve is the exact result. 



One the other hand the approximate kinetic description is very good in the purely 
irreversible case, 1/K = 0, (but it still misses some features, as discussed above and in 
Ref.Q). Then Eqs. © and (JUJ) simplify to 



$(t) (11 



dp 
dt 

5$ 1 _ e -(y(*)+*M) 

— = -2$(t) , (19) 

dt y(t) + z(t) 



where y and z are related to p and $ by Eqs. (JBJ) and 0. It is easy to show that at 
long times z(t) goes to a finite limit whereas y(t) goes to oo as t. As a result the kinetics 
approaches a non-trivial jamming limit with an algebraic 1/t behavior. The numerical 
solution of the approximate equation shows that the saturation density at the jamming 
limit pfl ~ 0.7422 is very close to the exact value, p^L = 0.74759.... The overall 
agreement with the exact density evolution is very good, as shown in Fig. 01 

5. Memory and Kovacs Effects 

Finally, we apply the approximate kinetic equations to the description of memory ef- 
fects. By introducing $ as the second state variable, one expects to obtain a response to a 
sudden change of the tapping strength that captures these memory effects. This is indeed 
illustrated in Fig. 0^ that shows how the density evolves when K is changed at t = 60 
from K = 500 to K = 5000 and from K = 500 to K = 200. At longer times, note that the 
density curve corresponding to the change from K = 500 to K — 200 crosses two times 
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Figure 4. (a) Approximate time evolution of the density for three cases: K = 500 (full 
curve), K = 500 from t = to t = 60 and K = 200 for t > 60 (dotted curve) and K = 500 
from t = to t = 60 and K = 5000 for t > 60 (dashed curve), (b) Parametric plot of 
insertion probability $ versus density for the three above cases. For completeness, the 
equilibrium insertion probability is added (dash-dot curve). 



the density curve calculated when K = 200 is kept constant along the process. This is the 
signature of a Kovacs-like effect which is more clearly displayed in the parametric plot of 
Fig. EJo: instead of going monotonically to the final value, the density decreases first and 
passes by a minimum after increasing again. This non-monotonic behavior of the density 
(or the volume) is precisely the Kovacs effect. The same phenomenon is illustrated for 
different values of the tapping strength in Figs and b: K is switched from 5000 to 
500 or from 5000 to 2000 after a time t = 250. In the former case, the density reached 
at t = 250 when K is changed is higher than the final (equilibrium) density reached at 
K = 500; nonetheless, the density starts by first decreasing to a value less than the final 
density until it reaches a minimum and finally increases again. On the other hand, no 
minimum is observed in the latter case. 

The Kovacs effect is traditionally represented as a "hump" in the volume as a function 
of time 27, 3(1 31]. The position and the height of the hump vary with the waiting time, 
i.e., the time spent at the initial tapping strength (recall that the process of compaction 
is always out of equilibrium), and with the amplitude of the shift in the tapping strength 
(or in other glassy systems, the shift in temperature). The same behavior is found here. 
It is illustrated in Fig. where we show the evolution of the inverse density for different 
protocols to reach the equilibrium steady state at K = 500: in the upper curve, K is 
switched from 5000 to 500 at t w = 240, in the intermediate curve K is switched from 2000 
to 500 at t w = 169 and in the lower curve K = 1000 at t w = 139. The waiting time t w is 
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Figure 5. (a)Approximate time evolution of the density for three cases:fT = 500 (lower 
curve), K = 5000 (dot-dashed curve), K = 5000 from t = to t = 250 and K = 500 for 
t > 250 (full curve) and K = 5000 from t = to t = 250 and K = 2000 for t > 60 (dotted 
curve). (b)Parametric plot of the insertion probability $ versus density for the four above 
cases: For completeness,the equilibrium insertion probability is added (dashed curve). 



chosen that the density reached at that time is equal to final density, p eq {K = 500). One 
observes that the height of the hump increases as the amplitude of the shift increases. 
The Kovacs effect, also observed in the simulation data (not shown here), results from 
competing trends that are qualitatively well accounted for by the present description 
in terms of two thermodynamic parameters. Including $ in the statistical-mechanical 
treatment is crucial for reproducing such behavior. 

6. Conclusion 

We have applied a two-parameter statistical-mechanical formalism inspired by the work 
of Edwards and coworkers to describe the compaction kinetics of the parking-lot model. 
The approximation gives a fair description of the kinetics, although it underestimates the 
relaxation time characteristic of the final approach to equilibrium. Inclusion of a second 
thermodynamic parameter allows one to qualitatively reproduce experimentally observed 
memory effects in granular compaction and to predict a Kovacs effect similar to what is 
observed in many glassy systems. 
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